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INTRODUCTION 


1. 1 Background 

For  the  past  few  years,  the  Air  Force  Rocket  Propulsion  Laboratory 
(AFRPL)  has  sponsored  a series  of  theoretical  and  experimental  programs  on 
the  retrieval  of  plume  flow-field  properties  by  analysis  of  the  infrared  radia- 
tive and  absorptive  properties  of  plumes.  This  report  covers  the  third  phase 
of  study  in  these  programs  by  The  Aerospace  Corporation.  The  first  phase1 
was  a study  of  the  classic  problem  of  retrieving  radial  profiles  of  gas  tempera- 
ture and  concentration  in  cylindrically-symmetric,  gaseous  plumes  from 
transverse  profiles  of  emission  and  absorption  (E/A  profiles)  obtained  in  a 
fixed  spectral  bandpass.  The  E/A  profiles  are  defined  in  terms  of  the  radial 
profiles  of  pressure,  temperature,  and  concentration  (pTc  profiles)  by 
integral  equations  of  radiative  transfer.  Retrieval  of  the  pTc  profiles  from 
the  E/A  profiles  involves  a numerical  inversion  of  these  integral  equations. 

In  this  study,  an  inversion  procedure  was  developed  and  incorporated  into  the 
computer  code  EMABIC.  This  code  was  used  extensively  in  the  present  study. 
The  inversion  algorithm  is  an  iterative  Abel  inversion.  The  well-known  Abel 
inversion  procedure  is  valid  for  optically  thin  sources;  for  the  general  case 
of  optical  thickness,  an  iterative  procedure  is  required.  The  code  has  been 
applied  to  several  synthetic  and  experimental  data  and  performs  satis- 
factorily as  a diagnostic  for  most  gas-only  plume  problems.  Some  problems 
occur  when  the  temperature  profile  has  a deep  minimum  on  the  plume  axis  or 
when  the  input  E/A  profiles  are  particularly  noisy,  even  if  they  are  adequately 
smoothed.  This  is  an  inherent  feature  of  inversion,  however,  and  is  not 


1.  S..  J.  Young,  Inversion  of  Plume  Radiance  and  Absorption  Data  for 
Temperature  and  Concentration,  AFRPL-TR-78-60,  U.  S.  Air  Force 
Rocket  Propulsion  Laboratory,  Edwards  Air  Force  Base,  Calif. , 

29  September  1978. 
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restricted  to  the  method  of  inversion.  A similar  inversion  code  has  been 

2 

developed  at  AEDC.  Recently,  a random  error  propagation  routine  was 

added  to  EMABIC  so  that  retrieval  error  could  be  estimated  automatically 

3 

from  E/A  measurement  error. 

The  second  phase  of  study  was  to  consider  multispectral  inversion 
and  the  effects  of  particle  loading  in  tactical  motor  plumes.  In  multispectral 
inversion,  retrieval  is  made  on  the  basis  of  how  E/A  spectra  vary  in  wave- 
length for  a fixed  measurement  line  of  sight.  It  was  found  that  this  inversion 
scheme  is  not  applicable  in  the  infrared  on  either  a monochromatic  or  wide 
band  spectral  scale  near  the  exit  plane  for  small  plumes  with  mild  temperature 
gradients  such  as  those  characteristic  of  tactical  rocket  motor.  Even  under 
the  most  ideal  circumstances,  temperature  and  concentration  retrieval  errors 
up  to  30%  were  encountered.  The  failure  of  the  method  is  caused  by  the  lack  of 
spatial  resolution  inherent  in  the  inversion  weighting  functions.  Results  of  this 
study  are  reported  in  Ref.  4.  Because  this  method  failed  for  purely  gaseous 
plumes,  it  was  never  applied  to  two-phase  plumes.  The  decision  was  made  to 
revert  to  the  multiposition  inversion  diagnostic  of  the  first  study  phase  and 
to  pursue  its  application  to  two-phase,  tactical  rocket  motor  plumes. 


2.  C.  C.  Limbaugh,  W,  T.  Bertrand,  E.  L.  Kiech  and  T.  G.  McRae, 

Nozzle  Exit  Plane  Radiation  Diagnostics  Measurements  of  the  Improved 
Transtage  Liquid  Rocket  Injector  Program,  AEDC-TR-79-29,  ARO  Inc. , 
Arnold  Engineering  Development  Center,  Arnold  Air  Force  Station,  Term. , 
March  1980. 

3.  S.  J.  Young,  Random  Error  Propagation  Analysis  in  the  Plume  Diagnostic 
Code  EMABIC,  AFRPL  Technical  Report,  to  be  published. 

4.  S.  J.  Young,  Multicolor  Inversion  Diagnostic  for  Tactical  Motor  Plumes, 
AFRPL-TR-80-30,  U. S.  Air  Force  Rocket  Propulsion  Laboratory, 

Edwards  Air  Force  Base,  Calif.,  May  1980. 
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1. 2 Scope  of  Present  Study 

The  concern  of  the  present  work  was  to  determine  how  applicable  the 
gas-only,  multiposition  inversion  diagnostic  code  EMABIC  is  to  plumes 
containing  particles.  The  work  focused  on  tactical  rocket  motors  where  the 
particle  loading  level  is  small,  that  is,  to  motors  where  the  particulate 
material  is  added  to  the  fuel  only  as  a stabilizer  (e.g. , Al^O^)#  or  to  motors 
where  the  plume  particulate  results  from  chemical  reactions  (e.g..  Carbon) 
but  not  to  motors  in  which  the  major  fuel  is  itself  a metal. 

The  general  procedure  of  the  study  was  to:  (1)  assume  reasonable  pTc 
profiles  and  particle  loading  levels  for  plumes  of  interest,  (2)  generate  E/A  j 

profiles  taking  account  of  particle  emission/absorption/scattering,  (3)  invert  \ 

these  E/A  profiles  with  the  inversion  code  EMABIC  as  though  they  were  for  a \ 

'purely  gaseous  plume,  and  (4)  compare  the  retrieved  profiles  with  the  known  j 

input  profiles.  The  results,  presented  later,  will  show  that  the  method  cannot  j 

tolerate  much  particulate  loading  without  introducing  appreciable  error  into  j 

1 

the  retrieved  gas  concentration  profile.  Temperature  retrieval  is  less  : 

sensitive  to  particle  loading.  A correction  procedure  was  derived  that  sub-  1 

stantially  increases  the  ability  to  retrieve  the  correct  profiles.  In  an  actual 
experiment,  E/A  data  would  be  obtained  in  a bandpass  appropriate  to  the  gas  \ 

species  of  interest.  This  data  would  reflect  both  gas  and  particle  effects.  j 

By  collecting  data  off  the  gas  absorption  band,  one  could  obtain  data  that 
primarily  reflects  just  the  particle  effects.  This  off-band  data  could  be  used 
to  "correct"  the  on-band  data  to  give  E/A  profiles  closer  to  the  gas-only 
profiles  needed  by  the  inversion  program.  In  this  work,  the  of band  data 
was  approximated  by  the  E/A  profiles  caused  by  the  particles  alone,  and  only 
the  first-order  correction  of  subtracting  the  off-band  radiance  from  the 
on-band  radiance,  and  dividing  the  on-band  transmittance  by  the  off-band 
transmittance  was  employed.  This  procedure  is  valid  only  if  the  plume  is 
optically  thin  to  both  gas  and  particle  absorption. 
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The  most  restricting  feature  of  the  present  work  was  the  use  of  the 
single-scattering  approximation  in  the  calculation  of  E/A  profiles.  The  use 
of  this  approximation  is  partly  justified  by  the  low  loading  levels  considered. 
These  calculations  are  considered  in  Section  2.  hi  Section  3.  the  optical 
properties  of  all  the  gases  (H^O.  CO^,  CO  and  HC1)  and  particles  (Al^O^, 
Carbon,  and  ZrOg)  considered  in  the  work  are  discussed,  and  in  Section  4 
specific  plume  models  are  considered. 
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2.  PLUME  RADIATION  CALCULATION 


2. 1 Single-Scattering  Band  Model 

In  this  section,  the  band  model  equations  £or  radiative  trans£er  in  a 
general,  nonuniform,  two-phase  (gas  and  particles)  medium  are  derived. 

The  formulation  considers  only  single  scattering,  but  allows  different 
temperatures  for  the  gas  and  particle  phases  and  treats  the  scattering  of 
external  source  radiation  (e.  g. , nozzle  exit  plane  radiation)  incident  on  the 
medium.  The  geometry  of  the  problem  is  shown  in  Fig,  1,  The  monochro- 
matic transfer  equation  for  the  radiance  along  the  primary  line  of  sight  (LOS) 
s directed  toward  the  sensor  is 


“ - = - Cor(s)+p(s)+  #(v,  s)]  N(v,  s)+a(s)  B (s)  +#(  v,  s)  B (s) 


4ir 


J P(s, 


(1) 


0)N{v,s,0)dO 


4tr 


where 


and 


m = wavenumber, 

O = shorthand  notation  for  direction  described  by  scattering 

angle  8 and  aximuthal  angle  0 , 

N(v,  s,0)  = radiance  directed  along  -0  direction, 

B .,(s)  = Planck  function  evaluated  at  gas  temperature  T (s), 

B (s)  = Planck  function  evaluated  at  particle  temperature  T (s), 
P P 

Of(s)  = particle  absorption  parameter, 

P(s)  = particle  scattering  parameter, 

p(s,0)  = scattering  phase  function  for  scattering  from  -0, 
direction  to  -s  direction, 

K(v,  s)  = gas  absorption  parameter. 
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The  parameters  or(s),  p(s),  and  p(s,0)  are  defined  by 


a{s)  ^ aa(s)  Np(s) 
p(s)  = CTs(s'}  Np(s) 

where  cra(s),  crg(s),  and  dcrg(0)/dn  are,  respectively,  the  absorption,  total 
scattering,  and  differential  scattering  cross  sections  for  the  particles 
(scattering  by  gas  molecules  is  neglected).  The  scattering  phase  function 
p(s,  0)  is  normalized  by 

ir 

4~  J p(B,0)d0  = ~J  p (8,d)  sine  de=  1. 

4ir  0 

Np(s)  is  the  particle  number  density.  The  gas  absorption  parameter  is 
defined  by 


#(v,  s)  = c(s)  p(s)  k{v,s) 

where  c(s)  is  the  gas  concentration  (mole  fraction),  p(s)  is  the  total  gas 
pressure,  and  k(v,s)  is  the  monochromatic  absorption  coefficient. 

Equation  (1)  is  written  for  a small  spectral  interval  Av  within  which  all 
parameters  are  independent  of  v except  for  K(v,  s)  and,  of  course,  N(v,  s) 
and  N(v,  s,  0). 

The  single-scattering  approximation  results  by  replacing  N(v,  s,  O) 
with  the  value  it  would  obtain  in  the  absence  of  scattering,  i.  e, , 
with  |3e0.  Specifically,  N(v,  s,fj)  is  approximated  as  the  solution  of 
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- dN(v^i,^.Q).  = . [a(a)  + K(v,a)]N<v,s,a,0) 


(2) 


+ <*(a)  Bp(cJ)  + k (v,  a)  Bg(a) 

evaluated  at  a = 0 and  with  the  imposed  boundary  condition  N(v,  8,  cr,  Q)  = 
Ne(v,s,Q)  at  cr=L(s,0).  The  result  is 

L(s,0) 

N(v,s,0)=  | Ca(a) Bp(a)  + KCv^BgfanT^ajT^crJda 
0 

+ NE(v,s,0)Ta[L(s,0)]T^  [v,L(s,0)] 


where 


Tw(a)=exp 


U 

I 


a(a')da' 


(4) 


t^(v,  a)  = exP 


#f(v,  ar'Jdcr'  | . 


(5) 


Substitution  of  Eq.  (3)  into  (1)  and  carrying  out  the  solution  of  the  first 
order  differential  equation  subject  to  the  boundary  condition  N(v,  s)  = 0 at 
s = L yields  the  following  solution  for  the  emergent  radiance  N(v)=N(\>,  s=0) 
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I IJggJW 


N(v)  = 


1 


(V,  8) 


+ T 


»(s)Bp(8)  Ta(s)Tp(8)  T^(v,s)  -Bg(s)  T^S)  Tp (s)  — 5^ 

(8)  Tp(s)  ^ Jp(s,0)  J ^ £ <*(<?) Bp(<r)  Ta(a)  t^(v,s  + ct) 

v,  s + a)  I 

3c J 


dT  ( 

- B„(c)  T Jo) 


g'  ' a 


4ir  l •'O 

dc 


+ Ne(v,s,0)  [L(s,fl)]  [v,^  + L(s,0)] 


or 


dO 


ds 


where 


Tp(s)  = exp 


3 

-I 


p(s')ds' 


(6) 


(7) 


and 


T^(v,8  + L)=  tk{v,s)  Tk(v,L). 


(3) 


In  this  form,  there  are  no  products  of  spectrally  correlated  functions,  and 
the  solution  for  the  mean  radiance 


H 


N = -57-  I N(v)  dv 
A V 


(9) 


can  be  obtained  immediately  by  replacing  all  references  to  t^(v,  s)  and 
di^  (v,  s)/dsinEq.  (6)byT(s)  and  dT(s)/ds,  respectively,  where  t (s)  is  an 
appropriate  band  model  transmittance  function.  Also,  we  assume  NE(v) 
independent  of  v and  replace  it  with  N„.  In  the  statistical  band  model,  the 

■k  1 

transmittance  derivative  can  be  written  as 


= - T(s)  k{b)  y(s) 


(10) 
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where  K(s)  is  the  mean  absorption  coefficient 


K(s) 


= mT  (*<'*- 

Av 


s)dv 


(11) 


and  y(s)  is  an  equivalent  width  derivative  function  whose  explicit  form 
depends  on  the  absorption  line  shape  and  the  approximation  used  to  account 
for  nonuniformities  along  the  optical  path.  In  effect,  this  function  accounts 
for  the  deviation  of  band  absorptance  from  Beer's  and  Lambert's  laws. 
Lamberts'  law  obtains  for  y = 1.  Use  of  Eq.  (10)  in  the  mean  radiance 
version  of  (6)  allows  the  final  solution  to  be  written  in  the  more  usable  form 


ds 


L 

r 

1 ^(8)^(8)  T(S) 

0 

Qt(s)  + Qg(s) 

Qt(s)  = or(s)Bp(s)  + K(s)  y(s)  Bg(8) 


Qo(s) 


= ^ fp(s, 

J4tt 


0) 


T [L(s,  Q)]  lLs+y^fl).3,N  (s,  o 

* T(S)  E 


L(s,r.) 

+ f T «r)  ^ 

J0  a T(S) 


<*(<*)  Bp(o)  + k(ct)  y(a  + s)  Bg(o) 


do 


(12) 


dO  . 


Q,p(s)  is  a thermal  source  term  and  Qg(s)  is  a scattering  source  term. 
The  transmittance  through  the  medium  along  the  primary  LOS  is 


FT  = VL)  VL)  T<L> 


(13) 
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Equations  (12)  and  (13)  are  the  general  £ormulas  used  to  compute  plume 
E/A  profiles.  Their  application  to  the  plume  problem  is  considered  in 
the  next  section. 

The  gas  band  model  formulations  used  to  compute  the  mean 
transmittance  ~ and  the  derivative  function  y are  described  in  detail 
in  Ref.  1 . Briefly,  the  formulation  employs  the  Malkmus  random  band 
model  with  provisions  for  Lorentz,  Doppler,  or  Voigt  line  profiles. 
Nonuniformities  in  the  plume  are  handled  with  either  the  Curtis-Godson  (CG) 
or  derivative  (DR)  approximations. 

A significant  feature  of  the  result  for  radiance,  Eq.  (12),  is  that  it 
correctly  accounts  for  the  nonmultiplicative  nature  of  band  model  trans- 
mittances.  This  account  is  made  in  three  places:  First,  radiation  that  is 
emitted  at  position  o on  the  scattering  LOS  and  will  eventually  be  scattered 
into  the  primary  LOS  is  not  simply  multiplied  by  t(o)  to  get  its  contribution 
to  the  scattering  source  function  at  s.  Rather,  its  contribution  is  correctly 
computed  by  multiplication  with  t(s  + o)/t(s)  where  t(s+g)  is  the  trans- 
mittance for  the  path  s=0  to  s=s(a=0)  and  then  on  to  o=o  treated  as  a single, 
coupled  LOS.  Note  that  for  monochromatic  radiation,  where  the  multiplica- 
tive property  does  hold,  this  procedure  would  yield  t(v,  s+ct)/t(v,  s)  e t(v,  ct) 
as  expected.  Second,  this  same  account  is  made  for  the  scattering  of  external 
source  radiation  into  the  primary  LOS.  The  third  account  is  the  use  of 
y{s+a)  instead  of  just  y(o)  in  the  expression  for  Qg(s). 

2.2  Calculation  Procedures 


Calculation  procedures  and  the  application  of  the  results  of  Section  2. 1 

to  plumes  are  described  here.  The  final  results  have  been  incorporated  into 

5 

a computer  code  which  is  described  in  a separate  publication. 


5.  S.  J.  Young,  User's  Manual  for  the  Plume  Signature  Code  EAFROF, 
AFRPL-TR-81-0&,  U.S.  Air  Force  Rocket  Propulsion  Laboratory, 
Edwards  Air  Force  Base,  Calif.,  January  1981. 


The  geometry  of  the  problem  is  shown  in  Fig.  2.  The  plume  is 
assumed  to  be  cylindrically  symmetric  and  axially  uniform.  The  obser- 
vation plane  is  located  a distance  Lj  downstream  of  the  nozzle  exit  plane  and 
the  plume  extends  a distance  L^  beyond  the  observation  plane.  In  gas-only 
radiation  prediction  and  inversion,  no  consideration  has  to  be  given  to  plume 
length  because  all  radiation  observed  originates  in  the  observation  plane. 

In  two-phase  plumes,  account  has  to  be  made  of  radiation  scattered  into  the 
observation  line  of  sight  from  other  parts  of  the  plume.  In  order  to  assign  a 
length  L,^  to  the  model  plume,  the  following  procedure  was  used.  First, 
realistic  flow-field  calculations  were  performed  (by  AFRPL)  for  plumes  of 
interest.  From  these  results,  the  decay  of  pressure,  temperature,  and  gas 
concentration  down  the  plume  axis  were  determined.  Second,  for  these  pTc 
variations  (or  approximations  to  them)  the  gas -only  radiance  at  the  observation 
station  directed  along  the  plume  axis  toward  the  exit  plane  was  computed. 
Finall/,  a similar  calculation  was  made  for  the  axially  uniform  model  (with 
pressure,  temperature,  and  gas  concentration  on  axis  equal  to  the  flow-field 
calculation  results  on  axis  at  the  observation  station)  with  L^  adjusted  until 
the  radiance  was  equal  to  the  previous  result.  Specific  values  for  Lj  and  L^, 
are  indicated  later  for  actual  cases  considered.  L^  was  generally  2 to  3 cm 
and  40  to  60  cm. 

The  s-axis  is  the  primary  LOS.  In  Fig.  2,  the  LOS  through  a full 
plume  diameter  is  shown,  but  as  the  LOS  is  scanned  across  the  lateral 
extent  of  the  plume,  it  becomes  increasingly  shorter.  The  a- axis  is  the 
scattering  LOS  and  is  described  by  the  value  of  s at  which  it  branches  off 
the  primary  LOS,  the  scattering  angle  0,  and  the  azimuthal  angle  0.  The 
length  of  the  scattering  LOS  is  L(s,  9,  0)  and  reflects  the  termination  of  the 
LOS  either  on  the  plume  broundary,  the  nozzle  exit  plane,  or  the  plume  end 
plane.  If  the  scattering  LOS  terminates  on  the  nozzle  exit  plane,  account  is 
made  of  thermal  motor  radiation  scattered  into  the  primary  LOS.  This 
radiation  is  modeled  as  arising  from  a disc  at  the  exit  plane  with  uniform 
temperature  Tg  and  emissivity  >5 . 
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Four  integrations  are  explicitly  indicated  in  Eq.  (12).  The  first  is 
over  the  line  of  sight  s.  In  the  original  form  of  EMABIC,  this  integration 
was  transformed  into  a radial  integration.  With  the  inclusion  of  scattering, 
it  is  more  convenient  to  retain  s as  the  integration  variable.  The  integration 
grid  for  s is  taken  as  the  set  of  points  defined  by  the  intersection  of  a given 
transverse  LOS  with  the  radial  zone  boundaries  (as  is  illustrated  in  Fig.  3). 
The  LOS  integration  is  performed  for  each  transverse  position  d of  the  grid. 
The  number  of  transverse  and  radial  zones  used  in  the  numerical  routines  is 
indicated  by  N.  Values  of  N = 7 and  10  were  used  in  this  study.  In  the 
example  of  Fig.  3,  N = 6.  The  radial  coordinate  index  is  i,  the  transverse 
coordinate  index  is  j,  and  the  primary  LOS  coordinate  index  is  k.  This  last 
index  runs  from  1 to  2(NH-j)+l.  The  midpoint  value  for  k is  km=N+2-j. 
The  radial  index  for  a given  value  of  k is 


N+2-k 

ksk 

2(j-l)-N+k 

k > k 

1 

The  radial  index  is  used  to  compute  the  actual  radial  coordinate  r=(i-l)  R/N 
for  each  k so  that  an  interpolation  on  tables  of  radial  functions  can  be  made 
to  get  the  function  value  appropriate  to  k.  The  coordinate  value  s for  each 
k is 


- 

1/2 

r 

1/2 

s =§- 

NZ  - (j-1)2 

± 

2.  2 
(i-ir-c-ir 

- 

. 

where  the  + sign  is  used  for  ksk^  and  the  - sign  for  k^k^.  A trapezoidal 
quadrature  routine  is  used  for  the  s-axis  integration. 


22 


At  each  grid  point  along  a given  transverse  LOS,  several  calculations 
for  the  integrand  of  the  s-integral  must  be  made.  The  calculations  for 
Ta(s)>  Tp(s),  t(s),  and  QT(s)  are  straightforward.  The  calculation  for  Q£(s) 
involves  three  more  integrations.  Two  of  these  integrations  are  over  the 
scattering  angle  0 and  azimuthal  angle  0.  In  Eq.  12,  these  integrations  are 
indicated  by  the  single  integral  over  fl . The  integration  ranges  are 
O£0£tr  and  O£0s:27r.  A Simpson  quadrature  is  used  for  0 near  zero  and 
trapezoidal  quadrature  is  used  for  large  0 and  for  the  0 integration.  The 
2ir  range  for  0 is  divided  into  equal  size  segments.  Most  of  the  work 
reported  here  was  done  with  N^  = 16.  The  0 integration  grid  is  selected  on 
the  basis  of  how  the  phase  function  p(0)  varies  with  0.  If  p(0)  peaks 
sharply  around  0 = 0,  the  density  of  grid  points  around  0~O  is  made  large. 

If  p(0)  is  relatively  flat,  a more  uniform  grid  is  used.  For  most  of  the 
work  performed  here,  the  following  grid  was  used: 


0.  = 0,5,  15,25,35,45,60,90,  120,  150  and  180°. 


(14) 


Further  details  on  the  selection  of  these  0 and  0 grids  are  given  in 
Section  4.  1.  1, 

The  fourth  integration  is  over  the  scattering  LOS  a , Before  this 
integration  is  performed,  an  analysis  of  the  geometry  is  carried  out  to 
determine  where  the  scattering  LOS  terminates,  and  how  long  it  is.  The 
relevant  geometry  for  the  analysis  is  shown  in  Fig.  4.  A scattering  LOS  is 
defined  by  che  distance  s along  the  primary  LOS  where  it  branches  off,  the 
transverse  distance  z,  and  the  angles  0 and  0.  The  cylindrical  coordinates 
for  a point  a aiong  the  scattering  LOS  are 
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NOZZLE  EXIT  PLANE 


y 


CYLINDER 

WALL 


Fig.  4.  Plume  Scattering  Geometry. 
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For  Os0Srr/2  or  3tt/2  5 0 £ 2ir,  the  scattering  LOS  is  directed  upward  and 
intersects  either  the  nozzle  exit  plane  or  the  cylinder  wall  above  the  scanning 
plane.  In  order  to  test  this,  Eq.  (16)  is  solved  for  a with  Z(a)  set  equal 
to  Lj.  This  value  for  o is  then  substituted  into  Eq.  (15)  and  solved  for  the 
value  rg  that  the  scattering  LOS  would  have  at  the  exit  plane.  If  rg  is  less 
than  R,  the  LOS  ends  on  the  exit  plane  and  the  length  of  the  LOS  is 


L = 


sinS  cos0 


If  rg  is  greater  than  R,  the  LOS  ends  on  the  cylinder  wall  and  the  length 
is  the  solution  of  Eq.  (15)  with  r(a)  replaced  with  R,  that  is 


L = 


2a 


2 1/2 

(b  + 4ac)  -b 
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whe  re 


a 

b 

c 


L-  . C . & 

= cos  0 + sin  G sin  0 


2(Xq  cosG  + y-Q  sinG  sin0) , 

_,2  2 2 

R -xo  -^o 


2 _ R2  2 
"R  "r0 


For  rr/2  £0£  3it/2,  the  scattering  LOS  is  directed  downward.  An  analysis 

similar  to  that  above  is  made  to  determine  whether  the  LOS  ends  on  the  end 

plane  or  the  cylinder  wall.  In  this  case,  Eq.  (16)  is  solved  for  a vdth 

Z(o)  =L,  and  used  to  compute  r . If  r <R,  the  LOS  ends  on  the  end  plane, 
fa  c c 

and 


L 


sinG  cos0 


If  rg>  R,  the  LOS  ends  on  the  cylinder  wall,  and  its  length  is 


L = 


2a 


2 1/2 
(b  + 4ac)  +b 


where  a,b,  and  c are  the  same  as  before. 

When  the  length  of  the  scattering  LOS  has  been  found,  it  is  divided  into  If 

s 

equal  sized  intervals,  and  die  N + 1 radii  of  the  interval  boundaries  computed 

from  Eq.  15.  These  radii  are  then  used  to  interpolate  on  tables  of  radial 

functions  to  get  the  integrands  of  the  cr  integral  at  the  N + 1 points.  A 

s 

trapezoidal  quadrature  routine  is  used  for  the  a integration.  All  of  the 
work  reported  here  used  Ng  = 10. 
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Although  the  formulation  procedure  discussed  so  far  does  not  require 
it,  a significant  simplication  has  been  made  in  all  of  the  work  reported  here  -- 
it  has  been  assumed  that  the  particle  size  distribution  and  index  of  refraction 
are  constants  throughout  the  plume.  In  fact,  this  simplification  is  not  very 
realistic  in  many  cases,  and  more  accurate  account  of  the  variability  of  these 
parameters  throughout  the  plume  is  planned.  Then  p(s,  A),  a(s),  P(s),  and  or(cr) 
are  independent  of  s or  a . 
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3.  GAS  AND  PARTICLE  OPTICAL  PARAMETERS 


All  of  the  gas  and  particle  data  used  in  the  analyses  of  Section  4 are 
considered  here.  Four  gas  species  and  three  particle  species  were  con- 
sidered, but  not  all  possible  gas/particle  pairings  were  analyzed. 

3.  1 Gas  3and  Model  Parameters 


The  four  gas  species  considered  are  H20,  C02,  CO,  and  HC1.  The 

spectral  bandpasses  and  resolution  used  for  each  species  are  listed  in 

Table  1.  The  H-O  bandpass  is  one  of  those  used  in  the  recent  AEDC  E/A 
^ 2 

measurements  program.  The  HC1  bandpass  is  simply  located  at  the 
center  of  the  emission  band.  The  C02  emission  band  is  located  on  the 
steeply  rising  blue  wing  of  the  CO^  emission  band.  This  is  the  only  spectral 
region  that  is  free  of  both  CO  and  H20  emission  contamination,  while  at  the 
same  time  is  not  so  strongly  absorbing  as  to  preclude  seeing  into  the  plume. 
The  CO  bandpass  is  located  just  beyond  the  red  wing  limit  of  the  CO-, 
emission  band  to  preclude  contamination  from  that  species.  Consideration 
on  the  spectral  purity  of  the  H20,  C02,  and  CO  bandpass  is  made  in  Ref.  4. 

The  HC1  and  CO  band  model  parameters  were  taken  from  the  NASA 

handbook;^*  the  H-O  and  CO,  parameters  are  those  of  the  plume  signature 
7-l(T  c — 

code  ATLES.  The  absorption  parameter  k and  line  density  parameter 


6.  C.  B.  Ludwig,  W.  Malkmus,  J.  E.  Reardon,  and  J.  A.  L.  Thompson, 
Handbook  of  Infrared  Radiation  from  Combustion  Gases,  eds.  R.  Goulard 
and  J.  A.  L.  Thompson,  NASA  SP-3080,  Marshall  Space  Flight  Center, 
Huntsville,  Ala.,  1973. 

7.  S.  J.  Young,  Band  Model  Parameters  for  the  2.  7 -um  Bands  of  H,Q  and 
CO,  in  the  100-3000°K  Temperature  Range,  TR-6076(6970)-4,  Tne 
Aerospace  Corporation,  El  Segundo,  Cam. , 31  July  1975. 

8.  S.  J.  Young,  Band  Model  Parameters  for  the  4,  3-um  Fundamental 
Band  of  CO,  in  the  100-3000C'K  Temperature  Range,  TR-007^> 

(6754-0$j-l,  The  Aerospace  Corporation,  El  Segundo,  Calif.,  19  ifeb.  1976. 

9.  S.  J.  Young,  "Evaluation  of  Nonisothermal  Band  Models  for  H-O,  " 

J.  Quant.  Spectrosc.  Radiat.  Transfer  JJ3,  29-45  (1977). 

S.  J.  Young,  Description  and  Use  of  the  Plume  Radiation  Code  ATLES, 
TR-0077(2753-64)-3,  The  Aerospace  Corporation,  El  Segundo,  Calif.'  13  May  1977. 
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10. 


T for  the  four  species  are  shown  in  Figs.  5-8.  Line  broadening  includes 
pressure  broadening  and  Doppler  broadening.  The  model  is  the  same  as 
that  used  in  the  code  EMABIC,  * and  the  relevant  parameters  are  given  in 
Table  2. 


3.2  Particle  Optical  Properties 


Many  problems  are  involved  in  modeling  the  absorption,  emission, 
and  scattering  properties  of  plume  particulates,  and,  at  present,  any 
method  of  modeling  is  open  to  criticism.  The  principal  problems  concern 
the  chemical  composition  (and  homogeneity  of  composition),  shape,  particle 
size  distribution,  and  spatial  distribution  of  plume  particles. 

The  three  species  of  concern  in  this  work  are  alumina  (A^C^), 

carbon  (C),and  zirconium  oxide  (ZrOz).  The  particles  were  assumed  to  be 

spherical  (or  compact  in  form  and  randomly  oriented)  and  homogeneous  so 

that  Mie  theory  could  be  used  to  model  their  scattering  properties.  The  index 

of  refraction  of  the  bulk  material  was  used  in  the  calculations  even  though 

some  doubt  exists  as  to  whether  or  not  this  is  the  proper  value  to  use  for 

11 

particles  whose  size  is  comparable  to  the  wavelength  of  interest.  Neglect 
of  this  effect  is  to  a large  part  justified  by  the  inaccuracy  with  which  we  know 
even  the  bulk  index.  The  index  data  used  in  the  present  work  are  considered 
in  Section  3.  2. 1. 

In  addition  to  the  index  of  refraction,  Mie  scattering  calculations 
require  the  size  or  size  distribution  of  the  particles.  Again,  little 
quantitative  data  exists  for  size  distributions  occurring  in  plumes.  The 
distributions  used  in  this  work  are  discussed  in  Section  3.2.2. 


11.  A.  B.  Pluchino,  S.  S.  Goldberg,  J.  M.  Dowling  and  C.  M.  Randall 
"Refractive  Index  Measurements  of  Single  Micron-Sized  Carbon 
Particles,"  Applied  Optics  19,  3370-3372  (1980). 


Model  Parameters  for  H-O.  (a)  Absorption  Parameter; 
ne  Density  Parameter. 


With  the  index  of  refraction  and  size  distribution  data,  standard 
Mie  theory  algorithms  were  used  to  generate  the  absorption  cross  section 
g , the  total  scattering  cross  section  a.,  and  the  differential  scattering 

cl  S 

cross  section  da  (0)/dO.  These  results  are  given  in  Section  3.  2.  3.  The 

absorption  parameter  a , the  total  scattering  parameter  (3,  and  the 

scattering  phase  function  p(fi)  required  by  the  single -scattering  band  model 

were  then  obtained  according  to  the  formulas  of  Section  2.  1.  The  total 

particle  density  N required  for  this  calculation  was  either  treated  as  a 
P 

parametric  variable,  or  computed  from  typical  mass  or  mole  fraction  loading 
levels.  These  levels  are  discussed  later  in  Section  4.  The  Appendix  treats 
the  conversion  from  one  measure  of  particle  loading  to  another. 


3.  2.  1 Refractive  Index 

The  application  example  of  Section  4.  1 required  the  index  of  refraction 

of  Al-,0,  at  2.  51  |xm  and  at  a characteristic  plume  temperature  of  ^ 1000K. 

^ i 12 
According  to  the  data  compilation  of  Whitson,  the  index  for  these  conditit  ns 

_ 7 

for  purr  Al^O^  is  n =”  1 . 7 and  K 6 x 1 0 (the  complex  index  is  m = n - i K). 

The  imaginary  component  of  this  index  (K)  is  much  too  small  to  explain  th< 

observed  optical  effects  of  A^O^  in  plumes.  A tentative  value  of 

m = 2.  0 - 0.  001 8i  is  used  here.  This  value  was  obtained  for  samples  of  Al-,0., 

13  C * 

actually  used  as  burn  stability  additives  in  propellant  mixtures.  The 
increased  value  of  K provides  an  improved  correlation  with  the  observed 
radiation  of  Al^O^  in  plumes,  but  is  still  judged  somewhat  low.  In  the  present 
analysis,  arbitrary  imaginary  values  of  0.01  and  0.05  were  also  used.  These 
increased  values  can  be  rationalized  by  noting  that  even  a slight  amount 
of  contamination  coating  a small  particle  can  dramatically  increase  its 


12.  M.  E.  Whitson,  Handbook  of  the  Infrared  Optical  Properties  of  Al-,0, . , 

Carbon,  M„0,  and  ZrO,,  TR-0075{5548}-2,  The  Aero3pace  Corpora1ux>n. 
ETslfuhTo  1975. 

13.  T.  D.  McCay,  U.  S.  Air  Force  Rocket  Propulsion  Laboratory, 

Edwards  Air  Force  Base,  Calif.,  private  communication,  1980. 
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emissivity  (by  increasing  the  effective  value  of  K). 

The  example  application  involving  carbon  (Section  4.2)  requires  the 
refractive  index  at  3.  39  pm  and  '-'2000K.  The  value  used  here  is  n = 2 and 
K = 1.0.  ' Owing  to  experimental  uncertainty  and  ignorance  of  the  type 

of  carbon,  the  imaginary  part  of  this  index  may  range  from  '**0.  4 to  1.5. 

The  optical  constants  for  Zr02  are  not  well  known  in  the  spectral 
region  of  interest  and  for  the  characteristic  temperatures  of  plumes.  The 
principle  data  available  are  those  of  Dowling  and  Randall*  ^ where  n and  K 
were  determined  between  4.55  and  25  pm  at  300  and  573K.  The  data 
between  4.  55  and  5.  00  pm  were  averaged  to  jbtain 

n ~ 1.45  K . 144  at  300K 

n ~ 1.39  K ~ .074  at  573K. 

By  linear  extrapolation  on  the  n-data,  the  values 

n » 1.2  at  1300K 

n *-  1.0  at  2000K 


14.  A.  B.  Pluchino,  "Emissivity  Spectra  of  Composite  Microscopic  Particles,1 
Applied  Optics,  Z0_,  531-533  (1981). 

15.  G.  N.  Freeman,  C.  B.  Lundwig,  W.  Malkmus,  and  R.  Reed,  De velopment 
and  Validation  of  Standardized  Infrared  Radiation  Model  (SIRRM)! 


were  obtained.  These  temperature  values  are  the  characteristic 
temperatures  of  the  plumes  for  two  cases  considered  in  Section  4.  3, 

A semilog  extrapolation  was  made  to  obtain 

K .01  at  1300K 

K - .0025  at  2000K 

These  n and  K were  used  for  all  the  spectral  bandpasses.  This  assumption 

is  partly  justified  by  the  fact  that  the  emissivity  and  absorptance  data  complied 
12 

by  Whitson  are  reasonably  constant  between  2 and  5 pm. 

A summary  of  the  refractive  index  data  used  in  this  work  is  given  in 
Table  3.  All  told,  the  quality  of  these  data  is  quite  crude  and  they  should  lot 
be  taken  as  definitive  values.  They  are,  however,  judged  adequate  to  display 
the  qualitative  features  desired  in  this  work.  The  index  for  ZrO for 
example,  represents  a particle  type  in  which  the  primary  optical  property  is 
scattering.  Little  absorption  or  emission  is  effected.  The  index  for  Al^C  3 
also  represents  a scattering  particle,  although  with  the  increased  values  for 
K,  it  .an  contribute  significantly  to  absorption  and  emission.  The  index  for 
carbon,  on  the  other  hand,  represents  a particle  type  in  which  absorption  md 
emission  are  at  least  as  important  as  scattering. 

3.  2.  2 Size  Distributions 


The  size  distribution  used  for  Al^O^  and  ZrC>2  *-s  a bimodal,  modified 
gamma  function  distribution  described  by 


f(r) 


alr 


+ 


c>h 


a2r 
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Table  3.  Particle  Index  of  Refraction, 


Particle 

n 

tc 

Reference 

a12°3 

2.51 

2000 

2.0 

0.0018 

0.01 

0.05 

12,  13 

Carbon 

3.39 

2000 

2.6 

1.0 

12 

Zr02 

2-5 

1300 

1.2 

0.01 

12,  16 

2000 

1.0 

0.0025 

40 


20.8 


4.  45  x 10 


= 35.92 


The  unit  of  r is  pm  and  the  unit  of  f(r)  is  pm"  . f(r)  is  normalized  by 
J f(r)  dr  = 1.  The  distribution  is  shown  in  Fig.  9.  The  distribution  is 

'0  17 

based  on  data  obtained  from  "small"  rocket  motors  and  is  characterized 
by  a large  submicron  peak  about  r = 0.  05  pm  and  a smaller  micron-size 
peak  about  r = 0.  5 pm.  Although  the  submicron  peak  is  an  order  of  magnitude 
larger  than  the  micron-size  peak,  the  scattering  parameters  are  comple’ely 
dominated  by  the  large -particle  distribution.  This  result  is  discussed  ir  the 
next  section. 

The  large  peak  of  this  distribution  is  reasonably  consistent  with  t e 

1 8 

mean  size  of  particles  deduced  by  McCay  et  al  from  visible  laser 

scattering  and  attenuation  measurements.  They  obtained  radii  values  of 
0.  30  to  0.  35  pm.  The  distribution  is  also  consistent  with  some  of  their  results 
for  ZrO In  two  of  their  tests,  they  obtained  mean  radii  of  ~ 0.  5pm.  In  a 
third  test,  they  obtained  results  which  could  be  interpreted  as  indicating 
particle  radii  of  0.9  to  1.5  pm.  In  order  to  simulate  this  latter  size,  a 
.nonodisperse  distribution  with  r = 1. 0 pm  was  also  used  in  the  present  work. 


17.  R.  Dawbarn,  Aluminum  Oxide  Produced  by  Solid  Rocket  Motors, 
Proc.  USAF/NASA  Int,  Spacecraft  Contamination  Coni!. , USAF 
Academy,  Colorado  Springs,  March  1978. 

18.  T.  D.  McCay,  W.  C.  Mundy,  D,  M.  Mann,  and  G.  S.  Meserve, 
Laser  Mie  Scattering  Measurements  of  Particle  Size  in  Solid  Rocket 
Motor  Exhausts,  JANNAF  12th  Plume  Technology  Meeting,  USAF 
Academy,  Colorado  Springs,  18-20  November  1980. 


The  particle  size  distribution  used  for  carbon  is  described  by  the 
modified  gamma  function 

a \ a 

f(r)  = a r e r 

with  a = l„  89x10**,  a = 4,  (3  = 32.57  and  y = 1/3.  The  unit  of  r is  |im 
and  of  f is  jim“*.  The  distribution  peaks  at  r = 0.  05  p.m  and  falls  to  a 
value  more  than  two  orders  of  magnitude  smaller  than  the  peak  value  at 
r = 0.5  pm.  The  function  is  shown  in  Fig.  10. 

3.  2.  3 Scattering  Cross  Sections 

19 

Standard  Mie  scattering  calculations  were  made  for  the  index  of 
refractions  and  size  distributions  just  discussed.  The  results  for 
are  shown  in  Table  4 and  Fig.  11.  These  results  are  the  sum  of  the  results 
for  each  mode  of  the  distribution  of  Fig.  9 taken  separately,  but,  for  all 
practicle  purposes, are  essentially  equal  to  the  results  for  the  large  size 
mode  alone.  The  small  size  mode,  although  roughly  two  orderp  of  magnitude 
larger  in  peak  number  density,  contributes  only  about  one  percent  to  the 
total.  The  total  cross  sections  are  tabulated  in  Table  4 for  the  three  assumed 
values  of  the  imaginary  part  of  the  refractive  index.  The  variation  of  K has 
little  effect  on  the  total  scattering  cross  section  og  and  almost  no  effeci  on 
the  extinction  cross  section  a = a + o . The  absorption  cross  section  , 

’3,  as  would  be  expected,  quite  sensitive  to  variation  in  K.  Except  possi  >ly 
for  the  largest  K , n is  much  larger  than  a and  indicates  that  particle 
scattering  should  be  more  important  than  particle  emission  and  absorption 
in  determining  the  role  of  particles  in  the  plume.  The  differential  scattering 
cross  section  dag(0)/df2  is  shown  in  Fig.  11.  In  the  important  forward 
diffraction  peak,  the  cross  section  is  relatively  insensitive  to  the  magnitude 
of  K . 

19.  D.  Deirmerdjian,  Electromagnetic  Scattering  on  Spherical  Polydispersions, 
Elsevier  Publishing  Co.,  Inc7,  New  York,  N.  Y. , 19^9. 


r*  - 11  * - - 


43 


Table  4.  Al^O^  Scattering  Cross  Sections. 


K 

°a(cm2) 

°s(cm2) 

oe(cm2) 

0.0018 

6.  39(-10) 

5.  86(-8) 

5.  92(-8) 

0.989 

0.  01 

3.  20(-9) 

5.  58{-8) 

5.  90(-8) 

0.946 

0.05 

1. l6(-8) 

4.  62(-8) 

5.  78(-8) 

0.  799 

The  results  for  carbon  are  shown  in  Fig.  12.  The  scattering  is  not 
as  sharply  peaked  in  the  forward  direction  as  it  is  for  Al-^O^  because  the 
mean  size  of  the  particles  (Fig.  10)  is  much  smaller.  The  result  a =**  cr 

S cl 

indicates  that  scattering  should  be  comparable  with  emission  and  absorption 
in  determining  radiation  transfer  in  carbon  bearing  plumes. 

The  total  scattering  cross  sections  derived  for  the  selected  n and  *' 
and  for  the  two  size  distributions  for  ZrC>2  are  tabulated  in  Table  5.  The 
differential  scattering  cross  sections  are  plotted  in  Figs.  13  and  14.  The 
principle  features  of  these  results  are:  1)  the  increasing  forward  scattering 
as  X decreases;  2)  the  strong  minima  between  100  and  130°  scattering 
angle  for  the  monodisperse  distribution  at  X = 2.51pm;  and  3)  most 
importantly,  the  extreme  smallness  of  cr^  and  dcrs/df2  for  the  high 
temperature  case.  This  latter  result  (Fig.  14)  is  a consequence  of  using 
n"“l  for  this  case.  Consequently,  at  2000K,  radiation  transport  is 
dominated  by  the  emission  and  absorption  properties  of  the  particulates 
and  is  essentially  independent  of  their  scattering  effects. 
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4.  EXAMPLE  APPLICATIONS 


The  analysis  procedure  outlined  in  Section  1.2  is  here  applied  to  the 
exhaust  plumes  of  three  model  tactical  motors.  The  first  exhaust  model  is 
representative  of  a minimum  smoke  solid  propellant  (MSP)  and  considers 
H^O/A^O^  as  the  gas/particle  phases.  The  second  is  representative  of 
an  advanced  liquid  propellant  (ALP  in  which  HCl/Carbon  is  the  gas/particle 
pair.  The  third  is  representative  of  a reduced  smoke  propellant  (RSP)  and 
considers  the  pairing  of  ZrC^  with  several  gas  species. 

4.  1 Minimum  Smoke  Propellant  (MSP)  (HoO/Al^O,)  Plume 

The  analysis  for  this  model  is  divided  into  two  sections.  In  the  first 
(Section  4.  1. 1),  a preliminary  MSP  model  is  considered,  and  analyses  are 
presented  on  numerical  accuracy,  relative  magnitudes  of  various  emisjion 
and  extinction  mechanisms  to  the  total  E/A  profiles,  and  compatibility 
between  the  codes  EAPROF  and  EMABIC.  In  the  second  (Section  4.  1. 2), 
the  final  analysis  for  this  example  is  presented. 


4.1.1  Preliminary  Analysis 

The  active  gas  is  and  the  particulate  is  Al^O^  . The  optical 

properties  of  these  species  and  the  spectral  band  (2,51  pm)  have  been 

discussed  in  Section  3.  For  all  of  these  preliminary  analyses,  the  particle 

scattering  data  for  K = 0.0018  were  used.  The  radial  pTc  data  for  the 

plume  model  are  given  in  Fig.  15.  The  gas  concentrations  are  constant  in 

radius.  The  particles  were  assumed  to  be  in  thermal  equilibrium  with  the 

gas  so  that  the  particle  temperature  profile  was  the  same  as  the  gas  tempera- 

ture  profile.  A flat  concentration  profile  of  Np  = 10  /cm  was  used.  The 

total  scattering  cross  section  (Table  4)  is  a = 5,  86  x 1 0"®  cm^  and  the  plume 

s 

diameter  is  D = 20  cm.  The  extinction  by  scattering  through  the  plume  is  then 
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1 - t = 1 - exp(-N  or  D) » 0. 1 1.  This  loading  level  is  thus  about  the  maximum 
p e 

that  can  be  tolerated  within  the  single-scattering  assumption. 

The  nozzle  exit  plane  was  modeled  as  a flat  disc  at  T = 1500K  and 
with  emissivity  e = 1.  The  scanning  plane  was  5cm  below  the  exit  plane, 
and  the  plume  extended  30cm  beyond  the  scanning  plane.  Except  for  the 
radial  temperature  profile,  the  plume  was  uniform.  The  Lorentz  line  profile 
and  Curtis-Godson  approximation  were  used.  The  number  of  radial  and 
transverse  zones  was  10.  Path  integrations  along  scattering  lines  of  sight 
were  done  with  the  path  divided  into  10  equal  length  segments  (regardless 
of  the  length  of  the  path). 

With  these  input  conditions,  numerical  experiments  were  carried  out 
with  the  code  EAPROF  to  determine  how  fine  an  angle  integration  grid  is 
required  to  obtain  accurate  scattering  results.  The  angle  variables  6 and  0 
are  defined  in  Fig.  4.  In  the  course  of  debugging  the  code,  a heuristic  grid 
of  1 1 points  was  used  for  the  scattering  angle  0.  (See  Eq.  14).  This  grid  is 
illustrated  in  Fig.  16.  The  function  plotted  is  the  weighting  function  sin0p(0) 
that  occurs  in  the  scattering  integral 
2ir  tt 

p(0)N(0,  0)  sin8d0d0  . 

0 0 

The  11  points  were  chosen  to  give  a good  coverage  of  this  function  with  a 
higher  concentration  of  points  in  the  region  where  the  function  is  largest. 

The  results  of  the  first  numerical  experiment  confirmed  the  suitability  of 
this  grid,  and,  consequently,  this  grid  is  used  throughout  the  remainder 
of  the  analyses.  In  this  analysis,  the  azimuthal  integration  grid  was  held 
at  N^  = 4(0  = O,  90°,  180°,  270°)  while  calculations  were  made  with  the 
following  scattering  angle  grids 
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0,  30,  60, . ..,  180  7 

0,  15,  30,„,.,  180  13 

0,  10,  20,.,.,  180  19 

0,  5,  10,...,  180  37 

The  results  are  shown  in  Fig.  17a„  Ng  is  the  scattering  component  of  the 

total  emission  at  the  transverse  position  z = 0.  That  is,  is  the  total 

radiance  minus  the  radiance  value  that  occurs  for  the  gas -only  case.  The 

dashed  lines  show  the  bounds  of  the  true  value  of  N , and  the  dot  shows  the 

s 

result  obtained  with  the  11 -point  scattering  angle  grid.  The  11  point  gives 
a value  of  N that  is  only  1.7%  higher  than  the  most  probable  true  value 
(the  center  of  the  dashed  band).  The  corresponding  value  of  total  radiance 
is  les3  than  1%  too  high. 

In  the  second  numerical  experiment,  the  scattering  angle  integration 

was  held  constant  at  N_  = 11  while  calculations  were  made  with  the  azimuth 

u 

angle  grids 

v’fe)  -3te) <N^(v) 

N0  = 4,  8,  16,  32. 

The  results  are  shown  in  Fig.  1 7b.  Again,  the  dashed  lines  indicate  bounds 

on  the  true  value  of  N . The  value  obtained  for  NM  = 16  is  4.  6%  too  low. 

s $ 

Combined  with  the  1.  7%  too-high  error  obtained  by  using  N.  = 11,  the  total 

_ 0 

error  estimated  for  N is  ~ 3%  too  low.  The  error  in  the  total  radiance  is 

s 

about  one  half  this  value.  This  order  of  magnitude  of  error  is  adequate  for 
this  work,  and  = 16  is  employed  in  the  remainder  of  the  work. 
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In  a second  analysis,  an  investigation  was  made  to  determine  the 
contributions  of  various  effects  to  the  total  radiance  and  extinctance  of  the 
plume.  The  results  are  shown  in  Tables  6 and  7 and  Figs.  18  and  19.  In 
the  tables,  z is  the  transverse  position  of  the  line  of  sight,  the  first 
column  is  the  radiance  (or,  in  Table  7,  the  extinctance)  for  the  gas-only 
case.  In  the  second  column,  the  particles  are  allowed  to  absorb  and  emit, 
but  not  scatter.  The  third  column  shows  the  full  results  for  absorbing, 
emitting  and  scattering  particles,  but  with  the  nozzle  exit  plane  emission 
turned  off.  Finally,  the  column  labeled  TOTAL  shows  the  result  obtained 
from  all  effects.  The  last  four  columns  show  the  contributions  (obtained 
by  various  substractions)  of  the  gas,  particle  emission/absorption,  particle 
scattering,  and  nozzle  radiation,  respectively. 

For  radiance,  the  principal  contributors  are  emission  from  the  gas, 
and  scattering  of  exit  plane  radiation  into  the  line  of  sight.  Emission/absorption 
by  the  particles  is  only  ~ 1%  of  that  due  to  gas,  and  scattering  of  radiation 
other  than  nozzle  radiation  is  even  smaller.  Note  that  the  scattering 
contribution  is  negative  for  z 4 cm,  Presumably,  the  net  effect  along  the 
longer  lines  of  sight  that  obtain  for  small  z is  that  more  radiation  is 
scattered  out  of  the  line  of  sight  by  particles  than  is  scattered  in. 

For  extinctance,  gas  is  again  the  major  attenuator.  Scattering  of 
radiation  out  of  the  line  of  sight  is  the  second  most  important  effect  and  is 
essentially  equal  to  the  gas  attenuation.  Absorption  by  particles  is  two 
orders  of  magnitude  down  from  the  gas  and  scattering  contributions. 

The  third  major  preliminary  analysis  was  a check  for  consistency 
between  the  plume  signature  code  EAPROF  and  the  diagnostic  code 
EMABIC  for  gas-only  plumes.  The  gas  and  calculation  conditions  previously 
defined  were  used  to  predict  gas  only  E/A  spectra  using  both  EAPROF  and 
EMABIC.  The  results  of  the  comparison  are  shown  in  Table  8.  The  percent 
difference  between  the  radiance  profiles  varies  from  ~0.  3 at  z = 0 to 
~ 13  at  z = 9 cm  (where  the  radiance  is  approaching  zero).  These  differences 
can  be  attributed  to  the  slightly  different  quadrature  routines  used  in  the  two 
codes.  The  differences  in  transmittance  are  even  smaller  (less  than  ~0.  L 

percent).  59 
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Fig.  19.  Plume  Extinctance  Contributions  for  the  Preliminary 
MSP  Plume  Model. 
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The  transverse  E/A  profiles  of  Table  1 were  then  inverted  using  the 
iterative  Abel  inversion  mode  of  EMABIC.  Figures  20  and  21  show  the 
convergence  characteristics  for  the  inversions.  Figure  20  is  the  rms 
difference  between  successively  retrieved  temperature  profiles  as  a function 
of  number  of  iterations.  Figure  21  shows  the  similar  result  for  concentration. 
The  significant  feature  of  these  curves  is  that  EMABIC  inverts  the  profiles 
generated  by  EAPROF  as  stably  as  it  inverts  its  own  transverse  profiles. 

The  accuracy  of  the  retrievals  is  shown  in  Ta.ble  9 and  Fig.  22.  In  Table  9, 
the  retrieved  temperature  and  concentration  profiles  from  the  20th  iteration 
are  compared  with  the  true  values.  The  differences  from  the  true  values  are 
plotted  in  Fig,  22.  Except  near  the  plume  boundary,  the  accuracy  of  retrieval 
for  either  the  EMABIC  or  EAPROF  profiles  is  quite  satisfactory.  The 
inversion  for  the  EMABIC  profiles  is  somewhat  better,  as  would  be  expected. 

These  results  justify  the  procedure  of  generating  E/A  profiles  with 
EAPROF  and  inverting  with  EMABIC.  We  can  be  assured  that  the  difference 
of  retrieved  profiles  from  true  profiles  for  cases  involving  particles  are 
indeed  caused  only  by  particle  effects. 


4.1.2  Final  Analysis 

The  final  analysis  for  the  MSP  plume  model  employed  the  pTc  profiles 
of  Fig.  23.  The  only  change  from  the  preliminary  analysis  is  in  the  tempera- 
ture profile  near  r = R.  The  nozzle  exit  plane  temperature  and  emissivity 
were  changed  from  the  preliminary  analysis  values  to  T = 800K  and  e = 0.  75. 


Particle  concentration  was  treated  as  a parameter  of  the  analysis  and  varied 


C *2 

from  zero  up  to  the  value  used  in  the  preliminary  analysis  (N  =10  /cm  ).  The 


scanning  plane  is  L^=  3 cm  from  the  exit  plane  and  the  plume  end  plane  is 
L,2  = 57cm  beyond  the  scanning  plane.  This  value  for  was  determined  as 
explained  in  Section  2.  2 with  an  axial  temperature  fall-oxx  of  dr  /dz  = -5K/cm. 
All  other  conditions  were  the  same  as  for  the  preliminary  analysis. 
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Table  9. 


Comparison  of  Retrieved  and  True  Radial  Profiles  for  the 
Preliminary  Minimum  Smoke  Propellant  Plume  Model. 


Temperature  Retrieval  (K) 

z(cm) 

True 

A 

EAPROF 

A 

0 

1500 

1500 

0 

1513 

13 

1 

1493 

1493 

0 

1505 

12 

2 

1470 

1471 

1 

1480 

10 

3 

1433 

1432 

-1 

1440 

7 

4 

1380 

1379 

-1 

1383 

3 

5 

1313 

1313 

0 

1314 

1 

6 

1230 

1228 

-2 

1225 

-5 

7 

1133 

1133 

0 

1127 

-6 

8 

1020 

1022 

2 

1012 

-8 

9 

893 

878 

-15 

860 

-33 

10 

750 

788 

38 

758 

8 

Concentration  Retrieval  (mole  fraction) 

z(cm) 

True 

EMABIC 

A 

EAPROF 

£ 

0 

0.  1500 

. 1500 

0 

. 1482 

0018 

1 

. 1500 

0 

. 1494 

0006 

2 

. 1498 

0002 

. 1501 

. 0001 

3 

. 1502 

. 0002 

. 1499 

00)1 

4 

. 1501 

. 0001 

. 1500 

0 

5 

1! 

. 1500 

. 0 

. 1496 

0034 

6 

1 

. 1507 

. 0007 

. 1503 

rO 

O 

o 

7 

. 1495 

0005 

. 1478 

00  12 

8 

. 1497 

-.  0003 

. 1486 

00' 4 

9 

. 1538 

. 0038 

. 1548 

. 0018 

10 

\ 

. 1465 

0035 

. 1451 

0049 

68 


(mole  fraction) 


-0.003 


-0.006 


\ / 
V 


i * 20 


■EMABIC 


EAPROF 


6\x 

\ ^ 


Convergence  Residuals  for  Inversion  of  >]/A  Profiles 
Generated  with  EMABIC  and  EAPROF  for  the 
Preliminary  MSP  Plume  Model. 
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For  these  conditions,  transverse  radiance /extinctance  profiles 
were  generated  with  K and  N as  parameters.  Example  profiles  for 

p 

K = 0.  01  are  shown  in  Figs.  24  and  25.  It  can  be  seen  that  the  effect 

of  particle  concentration  is  much  stronger  on  extinctance  than  on  radiance. 

The  variations  of  the  centerline  (z  = 0)  values  of  N and  t with  particle 

concentration  are  shown  in  Fig.  26  and  Table  10.  The  single  curve  for  1 - t 

results  for  all  three  values  of  k and  is  due  to  the  near  constancy  of  the 

total  extinction  cross  section  with  K . Radiance,  on  the  other  hand,  displays 

4 3 

a measurable  variation  with  k for  N greater  than  ~2xl0  /cm  . 

r 


A gas -only  inversion  of  the  transverse  profiles  was  made  with  the 

code  EMABIC.  Where  appropriate,  all  of  the  conditions  used  to  compute  the 

transverse  profiles  were  carried  over  into  the  inversion  (eg.  line  shape, 

number  of  zones,  etc. ).  Inversion  was  performed  in  the  iterative  Abel 

inversion  mode  with  no  presmoothing.  Convergence  was  deemed  complete 

when  the  rms  temperature  and  concentration  differences  between  successively 

iterated  profiles  were  less  than  10K  and  0.001,  respectively.  The  number 

of  iterations  required  to  meet  these  criteria  was  5 to  7.  Example  inversion 

results  for  K =0.  01  are  shown  in  Figs.  27  and  28.  With  increasing  particle 

density,  the  temperature  is  increasingly  underpredicted.  The  percentage 

errors  of  retrieval  at  the  plume  axis  (r  = 0)  are  shown  in  Table  11  and  Figs.  29 

3 3 

and  30.  Below  N ~10  /cm  , the  error  of  retrieval  is  negligible.  At  the 
P 

maximum  loading  level  consistent  with  the  single -scattering  approximation 
5 3 

(i.  e, , Np=  10  /cm  ),  the  temperature  is  underpredicted  by  15  to  20  percent 
(~250K)  and  the  retrieved  concentration  is  more  than  100  percent  too  high. 
These  results  are  nearly  independent  of  k within  the  range  of  variation  used. 
In  order  to  keep  the  concentration  error  within  acceptable  bounds  (say  10 
percent),  the  particle  loading  level  cannot  exceed  ~10  /cm  . 


The  first-order,  off-band  correction  procedure  described  in  Section  1.  2 

was  applied  to  the  cases  for  K =0.01.  The  gas-only,  par  tide -only,  and 

gas -plus -particle  radiance  and  transmittance  profiles  are  shown  in  Figs.  31 

4 3 

and  32,  respectively,  for  the  loading  value  Np  = 4xl0  /cm  . ! 
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Fig.  25.  Transverse'  Extinctance  Profiles  for  the  MSP  Plume  Model. 
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Table  11.  Temperature  and  Gas  Concentration  Retrieval  on  Plume 
Axis  for  Minimum  Smoke  Propellant  Plume  Model. 


N 

P-3 
(cm  ) 

K = 

0.  0018 

K ~ 

0.01 

K = 

0.  05 

T(°K) 

% Error 

T(°K) 

% Error 

T(°K) 

% Error 

0 

1506 

C 

I 

1506 

0 

1 

1506 

0 

100 

1506 

1506 

| 

1507 

j 

400 

1506 

1509 

1505 

1000 

1503 

1503 

\ 

1504 

4000 

1487 

-1 

1485 

-1 

1492 

-0. 

5 

10000 

1451 

-3 

\ 

1454 

-3 

1463 

-2. 

5 

40000 

1339 

-11 

1346 

-10 

1374 

-8 

100000 

1218 

-19 

1235 

-18 

1289 

-14 

N 

P 

(cm'3) 

K = 

0.  0018 

K = 0.  01 

K = 

0.  05 

C 

% Error 

c 

% Error 

C 

% Error 

0 

. 1501 

0 

1 

. 1501 

0 

. 1501 

0 

100 

. 1501 

. 1501 

. 1501 

400 

. 1504 

. 1494 

i 

. 1507 

1000 

. 1516 

. 1516 

1 

. 1517 

1 

4000 

.1558 

4 

. 1562 

A 

\ 

. 1567 

4 

10000 

. 1660 

11 

. 1663 

11 

. 1678 

12 

40000 

. 2175 

45 

. 2202 

47 

. 2253 

50 

1000000 

. 3306 

120 

. 3354 

124 

. 3473 

132 

78 


z (cm) 


. Transverse  Radiance  Profiles  for  Gas  and  Particle 
Phases  for  the  MSP  Plume  Model. 
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The  errors  in  approximating  gas -only  radiance  as  the  difference 

N **N  -N  are  shown  in  Fig.  33.  Figure  33a  shows  the  errors  as  a 
8 P 4 3 

function  of  z for  the  fixed  particle  density  of  4x10  /cm  while  Fig.  33b 

shows  the  error  at  z = 0 for  variable  particle  density.  Since  the  loading 

level  is  low,  the  error  of  the  approximation  is  not  large.  The  error  of 

approximating  the  gas-only  transmittance  by  t =t/t  is  zero  since  this 

& P 

relation  is  exact  in  the  single -scattering  approximation. 

The  accuracy  of  retrieval  using  the  off-band  correction  method  is 

shown  in  Fig.  34.  and  Table  12.  The  figure  shows  the  radial  error  profile 
4 3 

for  N =4x10  /cm  and  the  table  shows  the  centerline  error  as  a function 

P 5,3 

of  particle  density.  For  all  concentrations  up  to  10  /cm  loading,  the 

retrieval  error  of  temperature  and  concentration  was  less  than  the  con- 
vergence criteria  AT  = 10K,  Ac  = 0.  001.  Much  above  this  level,  it  is  expected 
that  even  the  off-band  correction  method  will  fail  because  both  the  gas  and 
particle  transmittance s through  the  plume  will  be  less  than  ~0.  9. 

Particle  loading  levels  for  various  retrieval  and  calculational  conditions 
are  given  in  Table  13.  The  conversions  between  the  three  measures  of  particle 
loading  are  described  in  the  Appendix. 


4.  2 Advanced  Liquid  Propellant  (ALP)  (HCI /Carbon)  Plume 


An  analysis  was  made  of  the  radiation  and  inversion  characteristics  of 
a model  two-phase  plume  containing  gaseous  HCI  and  carbon  particles,  The 
major  gas  species  of  the  plume  are  HCI,  HF,  and  N2  of  which  HCI  is 
considered  as  the  diagnostic  radiating  species  and  the  latter  two  as  back- 
ground gases.  The  plume  composition  is  given  in  Table  14.  The  gas 
temperature  profile  (Fig.  35)  is  parabolic  with  values  2278  and  1600K  at 
the  centerline  and  boundary,  respectively.  The  plume  radius  is  7cm.  The 
pressure  and  HCI  concentration  profiles  are  flat  with  values  p=  1.  18  at.m 
and  c = 0.  1 1 7. 


Table  12.  Temperat 


PARTICLE 


p = 1. 18  atm 
CHC1  ’ 117 


f 

1 


The  particle  temperature  profile  (Fig.  35)  is  parabolic  with  values 
2000  and  1600K  at  the  centerline  and  boundary,  respectively.  The  center- 
line  value  reflects  a 278  K lag  in  thermal  equilibration  between  the  particles 
and  the  gas.  Particle  concentration  was  assumed  to  be  flat.  The  nominal 
concentration  is  c^  = 0. 108,  but  this  parameter  was  varied  in  the  analysis. 

The  nozzle  exit  plane  was  modeled  as  a flat  disc  at  T = 1500  K and 
e = 0,  25.  The  scanning  plane  was  placed  2. 5 cm  below  the  exit  plane. 
Calculations  were  performed  with  the  Lorentz  line  profile  and  CG  approxi- 
mation. The  number  of  radial  and  transverse  zones  was  seven.  The 
N0=  11  and  = 16  grid  described  in  Section  4.  1. 1 was  employed.  Path 
integrations  along  scattering  lines  of  sight  were  done  with  the  path  divided 
into  ten  equal-length  segments. 

The  plume  was  assumed  to  be  axially  uniform  and  to  extend  83  cm 
beyond  the  scanning  plane.  With  this  length,  the  body  of  the  uniform  plume 
beyond  the  scanning  plane  yields  an  axially  directed  (toward  the  nozzle) 
gas -only  radiance  equal  to  that  emitted  by  a plume  having  a - 8.  57  K/cm 
linear  fall-off  in  temperature  along  the  axis. 

The  nominal  mole  fraction  loading  value  of  c =0.  108  converts 

P 8 3 

(see  Appendix)  to  a particle  density  loading  of  N^=  1.  37x  10  /cm  . The 

transverse  E/A  profiles  computed  for  this  loading  are  shown  in  Figs.  36 

and  37.  The  gas-only  and  particle-only  profiles  are  also  shown.  Clearly, 

for  this  loading,  the  plume  optical  properties  are  particle  dominated.  The 

effects  of  particles  are  ~17  times  those  of  the  gas  alone.  The  variations 

of  the  centerline  (z  = 0)  radiance  and  extinctance  with  particle  loading  are 

8 3 

shown  in  Fig.  38  and  Table  15  for  loading  up  to  4x10  /cm  . Most  of  the 
effect  of  the  parades  owes  to  their  emission  and  absorption,  but  a measurable 
amount  is  caused  by  scattering.  This  is  illustrated  in  Figs.  39  and  40  where 
the  radiance  and  extinctance  are  shown  as  computed  at  the  nominal  loading 
level  with  and  without  particle  scattering  (but  with  emission  and  absorption). 
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Table  15. 


Variation  of  Radiance  and  Extinctance  at  z = 0 
with  Particle  Density  for  Advanced  Liquid 
Propellant  Plume  Model. 


Np(cm-3) 

N(W/cm^-sr-ym) 

►—* 

1 

H 1 

0 

.09731 

.02666 

1(4) 

.09750 

.02674 

4(4) 

.09809 

.02703 

1(5) 

.09926 

.02758 

4(5) 

.1051 

.03035 

1(6) 

.1168 

.03587 

4(6) 

.1745 

.06297 

1(7) 

.2863 

.1149 

4(7) 

.7773 

.3345 

1(8) 

1.463 

.6238 

4(8) 

2.307 

.9783 
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Fig.  39.  Transverse  Radiance  Profiles  for  the  ALP  Plume  Model 
Computed  with  and  without  Scattering. 
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Inversion  of  the  total  (gas-plus)  E/A  profiles  of  Figs.  36  and 
37  yields  the  temperature  and  concentration  results  shown  in  Figs.  41 
and  42  (labeled  Retrieved).  The  retrieved  temperature  profile  is  300  to 
500  K too  low;  the  retrieved  HC1  concentration  is  two  orders  of  magnitude 
too  large.  The  variations  of  retrieval  error  at  r = 0 for  a range  of  loading 
levels  are  shown  in  Figs.  43  and  44  and  Table  16. 

The  off-band  correction  procedure  was  applied  using  the  particle-only 
profiles  shown  in  Figs.  36  and  37.  The  inversion  results  are  shown  in  Figs.  41 
and  42.  The  temperature  retrieval  has  been  worsened  (by  up  to  500  K at 
r = 0),  but  the  concentration  retrieval  has  been  improved  dramatically.  This  is 
in  keeping  with  previous  results  where  it  has  always  been  found  that  concentra- 
tion retrieval  error  is  much  larger  than  temperature  retrieval  error  and  is 
the  more  improved  with  the  application  of  the  off-band  procedure.  However, 
the  retrieved  concentration  is  still  a factor  of  two  different  from  the  true  value. 
The  conclusion  is  that  the  nominal  carbon  loading  of 


CP  = 

. 108 

11 

s* 

9.  32x  10  ^ g/cm 

N = 

1.  37x  108/cm3 

P 

is  too  large  for  application  of  the  gas -only  inversion  diagnostic,  even  with 
the  off-band  correction  procedure. 

Actually,  these  poor  results  are  to  be  expected  for  this  case.  In  the 
first  place,  the  nominal  loading  level  is  very  much  higher  than  the  level  at 
which  the  single -scattering  approximation  should  be  valid.  This  level  is 


°P 

. 0298 

m = 

2.  36x  10“°  g/cm 

P 

N = 

3.  47x  107  /cm3 

P 
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Fig,  42.  Concentration  Retrieval  Results  for  the  ALP  Plume  Model. 
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Temperature  and  Gas  Concentration  Retrieval  on  Plume 
Axis  for  Advanced  Liquid  Propellant  Plume  Model. 


'VO 


In  the  second  place,  the  first-order  off-band  correction  procedure  should 
not  necessarily  provide  much  improvement  (and  it  doesn't  for  temperature 
retrieval)  because  the  requirement  that  both  the  gas  and  particle  transmit  tances 
through  the  plume  be  > 0.  9 is  violated.  For  the  nominal  loading  level,  ti >e 
particle  transmittance  is  ^0.3  (Fig.  37).  The  maximum  loading  level  fo'-  which 
this  procedure  is  valid  is 


cp  = .0092 

m = 7.  55  x 10”^  g/cm^ 

p 

N = l.llxlO^/cm^ 

P 

4.  3 Reduced  Smoke  Propellant  (RSP)  (H^O,  CO.,,  CO,  HCl/ZrQ-,)  Plume 

Analyses  were  made  for  two  RSP  plumes  containing  H^O,  CO^,  CO, 
and  FC1  as  the  gas-phase  species  and  Zr02  as  the  solid-phase  species.  The 
gas  temperature  and  pressure  profiles  for  the  two  cases  are  shown  in  Figs.  45 
and  46.  Unlike  the  previous  plume  models,  temperature  and  pressure  here 
increase  with  radial  distance.  The  significant  differences  between  the  twc 
cases  is  the  lower  overall  temperature  for  the  model  of  Fig.  45.  This  model 
will  be  referred  to  here  as  the  "low-temperature"  model.  For  both  models, 
the  plume  radius  is  R = 5.  69  cm.  The  gas  concentration  profiles  were  assumed 
to  be  flat,  and  their  magnitudes  are  listed  in  Table  17.  Note  that  the  low  temper- 
ature model  does  not  include  HC1  as  an  active  species. 

Analyses  were  made  only  for  the  single  loading  level  of 
m -?.  28x10”  g/cm  for  each  model.  Two  spatial  distributions  were 

r 

considered.  The  first  is  simply  a uniform  distribution  (constant  in  r). 

The  second  is  the  Gaussian  distribution 


N(r)  r.  N e”Y^r^R^ 
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Fig.  46.  Radial  pTc  Profiles  for  the  High-Temperature  RSP 
Plume  Model. 
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This  distribution  is  shown  in  Fig.  47.  The  value  y = 9.  212  was  selected  to 
force  the  number  density  at  r = R/2  to  be  one-tenth  of  its  value  at  r = 0.  The 
centerline  value  is  related  to  the  mean  value  over  the  cross  section  of  the 
plume  by  NQ  = y N/(l-e~^). 

For  the  uniform  spatial  distribution,  the  number  densities  are 
(see  Appendix) 

3 

990/cm  (monodisperse) 

N = • 

P 3 

18000/cm  (bimodal) 

for  the  monodisperse  and  bimodal  size  distributions  (see  Section  3.  2.  2). 

For  the  Gaussian  distribution,  the  centerline  values  are 

3 

9120/cm  (monodisperse) 

1. 66  x 10^'.;m3  (bimodal) 

The  particle  temperature  profile  was  assumed  to  be  the  same  as  the 
gas  temperature  profile. 

The  nozzle  exit  plane  was  modeled  as  a flat  disc  at  T = 1460  and 
2080 K for  the  low-  and  nigh-temperature  cases,  respectively,  and 
e = 0.  25  for  both  cases.  These  temperature  values  are  the  gas  tempera- 
ture a at  r/R  = 0.  75.  The  scanning  plane  was  placed  3 cm  below  the  nozzle 
exit  plane.  Calculations  were  performed  for  the  Lorentz  line  profile  and 
CG  approximation.  The  number  of  radial  and  transverse  grids  was  10.  The 
Nq  = 11,  = 16  grid  described  in  Section  4. 1.  1 was  used.  Path  integrations 

along  scattering  lines  of  sight  were  done  with  the  path  divided  into  ten  equal 
length  segments. 
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The  plume  was  assumed  to  be  axially  uniform  and  to  extend  a distance 
L 2 beyond  the  scanning  plane.  The  values  employed  are  listed  in  Table  18. 
With  these  lengths,  the  body  of  the  uniform  plane  beyond  the  scanning  plane 
yields  an  axially  directed  gas -only  radiance  equal  to  that  emitted  by  a 
realistic  plume  having  a -5K/cm  linear  fall-off  along  the  axis  and  a linear 
pressure  increase  that  restores  the  plume  pressure  to  1 atm  at  the  point 
the  temperature  reaches  300  K. 

The  gas-only  E/A  profiles  are  shown  in  Figs.  48-51.  The  particle 

loading  is  so  small  that  essentially  these  same  profiles  (with  one  exception) 

result  for  all  combinations  of  spatial  and  size  distributions.  The  centerline 

(z  = 0)  results  for  all  conditions  are  given  in  Table  19.  The  one  exception 

occurs  for  the  bimodal  size  distribution  and  Gaussian  spatial  distribution  for 

the  H^O  bandpass  in  the  low-temperature  model.  The  gas-only,  particle-only 

and  total  radiance /extinctance  profiles  for  this  case  are  shown  in  Figs.  52 

and  53.  The  particle  density  and  mole  percent  loading  for  this  case  are 
4 3 

N =1.8x10  /cm  and  c =.  0022%,  These  are  mean  values  over  the  plume. 

P P 

The  results  of  Fig.  52  and  53  demonstrate  that  this  particle  loading 
influences  the  gas-only  E/A  profiles  only  slightly  and  only  near  z = 0. 

Retrieval  analysis  was  made  only  for  this  one  case.  Inversion  was 
performed  on  both  the  total  radiance /extinctance  profiles  of  Figs.  52  and 
53  as  well  as  the  profiles  resulting  from  the  first-order,  off-band  correction 
procedure.  The  results  are  presented  in  Tables  20  and  21  and  Fig.  54. 

The  maximum  temperature  retrieval  error  is  only  <~7%  even  without  the 
correction  procedure.  With  correction,  the  maximum  error  is  reduced  to 
<~3%.  Without  correction,  the  maximum  H^O  concentration  retrieval  error 
is  ~60%.  With  correction,  the  error  is  reduced  to  ~8%. 

Since  this  retrieval  case  represents  the  worst  case,  the  gas -only 
inversion  diagnostic  (with  first-order,  off-band  correction)  is  applicable 
to  all  cases. 
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Table  18. 


Effective  Plume  Lengths  for  Reduced 
Smoke  Propellant  Plxime  Models. 


Case 


Species  L(cm) 


Low- Temperature 
Model 


CO 

57 

co2 

74 

h20 

44 

High- Temperature 


rn 


o o.  2 0. 4 0. 6 0. 8 1.0 

z/R 


Fig.  48.  Transverse  E /A  Profiles  for  CO.  L = Low- Temperature 
Case,  H = High-Temperature  Case. 


Ill 


i 


Transverse  E/A  Profiles  for  HC1.  L = Low- 
Temperature  Case,  H = High-Temperature 
Case. 
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Transverse  Extinctance  Profiles  for  H^O/ZrO^  in 
the  Low-Temperature  RSP  Plume  Model  with 
Bimodal  Size  and  Gaussian  Spatial  Distributions. 
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Fig.  54.  Temperature  and  Concentration  Retrieval  Results 
for  the  Low-Temperature  H-O/ZrO-  RSP  Plume 
Model  with  Bimodal  Size  ana  Gaussian  Spatial 
Distributions. 


5.  summary  and  conclusions 


The  primary  goal  of  this  work  was  to  determine  quantitative  limits 
of  particle  loading  in  realistic  tactical  rocket  motor  plumes  for  which  the 
gas  properties  of  the  plume  could  be  retrieved  with  standard  E/A  inversion 
diagnostics  without  having  to  account  for  the  radiation  and  absorption  effects 
of  the  particles.  The  procedure  was  to  assume  flow-field  properties  for 
two-phace  plumes  of  interest,  generate  E/A  profiles  with  account  of  particles 
using  a single-scattering  plume  radiation  model,  retrieve  the  gas  properties 
from  the  E/A  profiles  with  the  gas -only  inversion  code  EMABIC  under  the 
assumption  that  the  profiles  were  caused  by  gas  alone,  and  compare  the 
retrieved  gas  properties  with  the  assumed  properties.  Generally,  the  degree 
of  particle  loading  was  treated  as  a parameter. 

Two  important  results  were  obtained.  The  first  is  vhat  the  limit  of 
particle  loading  at  which  reasonable  (<10%  error)  retrieval  results  can  be 
obtained  is  generally  smaller  than  the  nominal  loading  level  for  the  plume. 

The  important  implication  of  this  result  is  that,  to  the  extent  that  the  systems 
studied  here  are  typical,  E/A  diagnostics  on  plumes  generated  by  even  low- 
and  reduced-smoke  type  propellants  requires  some  account  of  particle  effects. 
The  second  result  is  that  the  maximum  loading  level  for  acceptable  gas 
temperature  retrieval  is  much  higher  (about  an  order  of  magnitude)  than  that 
for  gas  concentration  retrieval.  Consequently,  in  applications  where  tempera- 
ture retrieval  is  of  primary  concern,  the  use  of  gas-only  E/A  diagnostics 
may  be  justified  even  though  the  total  retrieval  results  may  be  substantially 
in  error,  These  loading  level  results  are  summarized  in  Table  22  for  the 
three  systems  studied  here. 

Analysis  was  also  made  of  a correction  procedure  in  which  first-order 
corrections  were  made  to  the  total  gas/particle  E/A  profiles  by  using 
particle -only  E/A  profiles  obtained  outside  the  gas  absorption  band.  The 
corrected  profiles  provide  better  estimates  to  the  gas-only  profiles  needed 


il  1 


\ 


in  the  gas -only  inversion.  The  particle  loading  limit  for  valid  use  of  this 
procedure  is  the  value  for  which  the  total  extinction  of  radiation  by  particles 
over  a full  diameter  of  the  plume  is  ~10  percent.  (Note  that  if  this  condition 
is  met.  then  the  condition  is  also  met  that  the  attenuation  by  scattering  alone 
over  this  path  is  less  than  ~10  percent.  The  latter  condition  is  required  by 
the  single -scattering  assumption  used  in  this  work. ) For  the  two  cases 
(MSP  and  RSP)  where  the  nominal  loading  limit  roughly  satisfies  this  condition 
(see  Table  22),  the  use  of  the  procedure  gave  retrieved  gas  properties  that 
were  accurate  to  within  the  convergence  criteria  set  on  the  inversion.  For 
the  ALP  model,  the  nominal  loading  level  was  well  above  this  limit,  and  the 
retrieval  results  were  poor. 
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Table  22.  Summary  of  Critical  Loading  Levels. 


Model 

Species 

_3 

Loading  Levels  (cm  ) 

Te>0.  9 

ec<10% 

eT<10% 

Nominal 

MSP 

ai2o3/h2o 

8.9(4) 

1.0(4) 

6.  0 (4) 

1.0(5) 

ALP 

Carbon/HCl 

1.1(7) 

1. 0(5) 

1. 0(6) 

1.4(8) 

RSP 

Zr02/H20(2) 

9.4(5) 

3.0(3){3) 

2.0(4) 

1. 8(4)(4) 

(1) 

(2) 

(3) 

(4) 


K=0.  01 

Bimodal  size  distribution,  Gaussian  spatial  distribution, 
low-temperature  model. 

4 3 

Linear  extrapolation  from  60%  error  at  = 1.8x10  / cm  . 
Average  over  observation  plane. 
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APPENDIX 


CONVERSION  BETWEEN  NUMBER  DENSITY,  MASS 
LOADING,  AND  MOLE  FRACTION 


Several  measures  of  particle  loading  in  plumes  are  employed.  In  this 

appendix,  the  formulas  are  derived  for  conversions  between  loadings  expressed 

-3  3 

as  particle  number  density  (cm  ),  mass  loading  (g/cm  ),  and  mole  fraction 
or  mole  percent. 


The  conversion  between  number  density  and  mass  loading  is  the 
simplest  and  is  independent  of  the  plume  gas  properties.  Let  Np  be  the 
particle  number  density,  V the  mean  volume  of  particles,  and  p the  material 
density.  Then,  the  mass  loading  is 

mp  = Np  V p . (Al) 


For  a monodisperse  particle  size  distribution,  the  mean  volume  is  simply 
V=4Tra  /3  where  a is  the  particle  radius.  For  a modified  gamma  distribution 
described  by 

1(a)  = Aa'V3aV, 


the  mean  volume  is 

y tt  f(a)  d a 

0 


4tt  A 

3 Yp  (a  + 4)/y 


r 


(A2) 
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where  r is  the  gamma  function.  If  the  loading  is  uniform  over  the  cross- 
sectional  area  of  the  plume  at  a fixed  axial  station,  Eq.  (Al)  can  be  applied 
directly  to  obtain  the  conversion  formula.  For  the  MSP,  ALIJ  and  one  case 
of  the  RSP  plume  models,  the  loading  was  assumed  to  be  uniform,  and  the 
resulting  formulas  are  given  in  Table  Al.  The  density,  mass,  and  volume 
data  required  for  the  derivation  are  given  in  Table  A2. 

The  RSP  plume  model  was  also  analyzed  with  the  Gaussian  spatial 
distribution 


N(r)  = Noe"v^R) 
Y = 9.212 


(A3) 


In  Eq.  (Al),  we  assume  that  a given  is  the  mean  value  over  a cross- 
sectional  area  and  replace  on  the  right-hand  side  by  its  mean  value  over 
the  area.  In  this  case 


N = fr.-V(WR) 


R 


I 


r e 


dr  = N 


l-e"Y 


(A4) 


The  conversion  formula  is  then 

mp 


(A5) 


This  result  for  the  relationship  between  and  Nq  for  the  RSP  plume 
model  with  the  Gaussian  spatial  distribution  is  given  in  Table  Al, 
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Table  Al. 


Particle  Loading  Conversion  Formulas.  Units  are 

-3  3 

cm  for  N or  N . g/cm  for  m . and  mole 
P o g 

fraction  for  c . 


MSP  Plume  Model 

(ai2o3/h2o) 

c = 1 

P .3 

1+1.16x10 

N = 1.  17  x 1012m 

P P 

m 

P 

ALP  Plume  Model 

(C/HC1) 

c = 1 

p -5 

1 + 7.69  x 10 

N = 1.47  x I013m 

P P 

m 

P 

RSP  Plume  Model  (ZrC>2/H20) 


1 + 1-  04  x 10 
m 


1 + 8.49  x 10 
m 


4.  34  x 1010  mp 
7. 82  x 1011  m 

P 

4. 00  xlO11  m 

P 

7.  27  x 1012  m 


low-temperature 

high- temperature 

monodisperse  size/uniform  spatial 
bimodal  size /uniform  spatial 
monodisperse  size/Gaussian  spatial 
bimodal  size /Gaussian  spatial 


% 


Table  A2. 


Density,  Mass,  and  Volume  Parameters  for 
Solid- Phase  Species. 


Species 

p(g/cm3) 

w(amu) 

?(pm3) 

A12°3 

3.  70 

102 

0.  230 

C 

2.  25 

12 

0.  030 

ZrOz 

5.5 

123 

0.  230  bimodal 

4.  19  monodisperse 

The  conversion  from  either  particle  density  or  mass  loading  to  mole 
fraction  requires  the  gas  properties  of  the  flew  field.  The  number  of  gas 
molecules  per  unit  volume  for  a gas  at  temperature  T and  pressure  p is 

ng  = Cp/T  (A6) 

21  3 

where  C=  7.34x10  molecules  K/atm-cm  . The  number  of  particle 
molecules  per  unit  volume  is 

N V p m 

nP  = -£r-  = dr  (A?) 

where  p is  the  particle  mass  density,  w is  the  atomic  weight  of  the  molecules 
comprising  the  particle,  and  f = 1. 67x  10-24 g/amu.  The  particle  mole 
fraction  loading  is  thus 


c 


P 


n 


n +n 
P g 


1 

1 +nwf 

-JS 

m 

P 


(AS) 


The  measure  of  mole  fraction  loading  used  in  this  report  refers  to 
mean  values  over  a cross-sectional  area  of  the  plume  at  a fixed  axial  station. 
The  modification  required  to  account  for  this  is  the  replacement  of  n as 
given  by  Eq.  (A6)  by  the  average  value 

"g  = 0 (A9> 

0 
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m 


1 


11 


l»' 


••.  t 

>.i 


fc- 


*'■ 


i 


b 


fe  * 

%j 


M 


For  the  MSP  and  ALP  plume  models  considered  in  Section  4,  pressure  is 
constant  in  r and  T(r)  is  reasonably  well  approximated  as  the  parabola 


T(r)  = T -(T  -T,,)  (r/R)‘ 


where  Tq  is  the  plume  centerline  temperature  and  is  the  plume 
boundary  temperature.  Equation  (A9)  can  then  be  integrated  directly  to 
give 


( 


r d r 


= Se 

0 T0-(T0-XB)(r/R)2  T® 


(A10) 


where 


T is  the  effective  temperature 


Te  = 


To-TB 


(All) 


For  the  MSP  model,  T =1500K,  T- “•  750  K and  T = 1082K.  For  the  ALP 

o x)  e 

model,  Tq  = 2278K,  Tg  = l600K  and  Tg  = 1919K,  The  conversion  formula 


resulting  from  these  values  of  Tg  in  Eq.  (A10)  and  the  resulting  n in  Eq.  (A8) 


for  the  MSP  and  ALP  models  are  given  in  Table  Al.  The  same  formulation 
was  used  to  obtain  the  conversion  formulas  for  the  two  RSP  plume  models 


except  that  n was  computed  by  numerical  integration  of  Eq.  (A9)  since 
8 


neither  was  p(r)  constant  nor  T(r)  representable  as  a parabola.  The 
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